Next Article in Journal
Digital Quantum Simulation of Linear and Nonlinear Optical Elements
Previous Article in Journal
Information Geometric Perspective on Off-Resonance Effects in Driven Two-Level Quantum Systems
Previous Article in Special Issue
Information Entropy Approach for a Disorderless One-Dimensional Lattice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shannon Entropy in Confined He-Like Ions within a Density Functional Formalism

Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Mohanpur-741246, Nadia, WB, India
*
Author to whom correspondence should be addressed.
Quantum Rep. 2020, 2(1), 189-207; https://doi.org/10.3390/quantum2010012
Submission received: 6 January 2020 / Revised: 3 February 2020 / Accepted: 4 February 2020 / Published: 10 February 2020
(This article belongs to the Special Issue Applications of Information Entropies in Quantum Science)

Abstract

:
Shannon entropy in position ( S r ) and momentum ( S p ) spaces, along with their sum ( S t ) are presented for unit-normalized densities of He, Li + and Be 2 + ions, spatially confined at the center of an impenetrable spherical enclosure defined by a radius r c . Both ground, as well as some selected low-lying singly excited states, viz., 1sns (n = 2–4) 3S, 1snp (n = 2–3) 3P, 1s3d 3D, are considered within a density functional methodology that makes use of a work function-based exchange potential along with two correlation potentials (local Wigner-type parametrized functional, as well as the more involved non-linear gradient- and Laplacian-dependent Lee-Yang-Parr functional). The radial Kohn-Sham (KS) equation is solved using an optimal spatial discretization scheme via the generalized pseudospectral (GPS) method. A detailed systematic analysis of the confined system (relative to the corresponding free system) is performed for these quantities with respect to r c in tabular and graphical forms, with and without electron correlation. Due to compression, the pattern of entropy in the aforementioned states becomes characterized by various crossovers at intermediate and lower r c regions. The impact of electron correlation is more pronounced in the weaker confinement limit and appears to decay with the rise in confinement strength. The exchange-only results are quite good to provide a decent qualitative discussion. The lower bounds provided by the entropic uncertainty relation hold well in all cases. Several other new interesting features are observed.

1. Introduction

A particle in an impenetrable box of infinite height has served the role of a simple, elegant pedagogical tool to illustrate the effects of boundary condition on the energy spectrum of a quantum system. Understanding of such a system in some sub-region Ω of space (in contrast to the “whole” space available in a free system) offers new insights to simulate realistic situations in highly inhomogeneous media or in an external field. Matter constricted under such an extreme pressure environment gives rise to a wide range of novel changes (from the respective free counterpart) in energy spectra, electronic structure, chemical reactivity, ionization potential, polarizability etc., depending on the geometrical forms of the cavity and dimensions. This has inspired a variety of theoretical and experimental works. Some prominent applications are found in the context of the cell model of a liquid, superlattice structure, quantum dot, quantum wire, atoms/molecules encapsulated inside nanocavities (like fullerene, zeolite sieves, porous silicon, carbon nanotube), modeling defects in solids, confined phonons (or plasmons, polaritons, gas of bosons), as well as astrophysical phenomena such as the mass-radius relation of white dwarfs, ionized plasma etc. The topic is vast, and there has been a burgeoning growth of activity as evident from an extensive literature having many excellent comprehensive reviews. The interested reader may refer to the following reviews [1,2,3,4,5] and the references therein.
The first report of a confined hydrogen atom (CHA) within a sphere having rigid impenetrable walls was published as early as in 1937 [6], imposing the Dirichlet boundary condition that the wave function vanishes at the boundary. Subsequently, many attempts have been made to estimate the eigenvalues and eigenfunctions invoking a wide range of approximate analytic, semi-analytic, and purely numerical schemes. Here we mention a few like Rayleigh-Schrödinger perturbation theory, the Wentzel-Kramers-Brillouin method, the power-series solution, hypervirial theorem, Padé approximation, Lie algebraic treatment, super-symmetric quantum mechanics, the Lagrange-mesh method, the asymptotic iteration method, searching the zeros of the hypergeometric function, the generalized pseudospectral (GPS) method and the Hartree-Fock (HF) method [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. In recent years, the exact solution of the Schrödinger equation has been found in terms of the Kummer-M function (confluent hypergeometric) [16]. As the boundary approaches the nucleus and the volume of confinement squeezes, one notices a monotonic increase in energy in CHA. Another interesting feature is that, in contrast to a free atom, due to the breaking of the symmetry in CHA, it is characterized by different energy eigenvalues, eigenfunctions, and reduced degeneracies. On the other hand, new degeneracies, namely simultaneous, incidental and inter-dimensional degeneracy, which are non-existent in the free system, are introduced in CHA afresh. Apart from the effect of compression on the ground and various energy levels, properties such as the dipole shielding factor, nuclear magnetic screening constant, hyper-fine splitting constant, pressure, static and dynamic polarizability, etc., were examined.
An analogous study of the compressed He atom is a prototypical non-trivial mathematical problem. Due to the presence of inter-electronic repulsion, the S U ( 3 ) symmetry of the simplified one-electron case is broken, which promises many exciting physics. Ever since the variational calculation of energy variation [23] with respect to the cage radius and function of pressure, vigorous attempts have been known. Some of them include the Roothaan-HF-type calculation with the Slater-type basis [24] or with its modifications [25], configuration interaction [26,27], quantum Monte Carlo [28], a host of direct variational schemes with the appropriate choice of cut-off function [8,29,30,31], variational method with B-splines basis [32], and Rayleigh-Schrödinger perturbation theory [30,33]. Some other prominent works are explicitly correlated Hylleraas-type wave functions within the variational framework [34,35,36,37,38,39,40,41], a combination of the quantum genetic algorithm and HF method [42], variational Monte Carlo [43,44], HF calculation employing local and global basis sets [45], and so on. Whereas a vast majority of publications have focused on the lowest state, low-lying excited states were also treated quite decently. For example, 1sns 1,3S states in [29,30,35,38,39,40,41,42,44,46], 1s2p 3P, 1P states in [29,42,46], singly excited 1s3d 3D, 1D and some doubly excited states in [42,46] etc., using a host of theoretical approaches giving results of varied accuracy.
All the above works pertain to the wave function-based methods. In the past two decades, some results have been published within the alternative density-based concept, the so-called density functional theory (DFT) [47,48,49]. Thus, within an exchange-only framework (using two exchange functionals, viz. local density approximation (LDA) [47] and Becke-88 exchange potential [50]), the desired Kohn-Sham (KS) equation was solved satisfying the Dirichlet boundary condition for many-electron systems via the numerical shooting method [51]. The usefulness of a one-parameter hybrid exchange functional (including a fraction of exact exchange and the Perdew-Burke-Ernzerhof functional) for treatment of confined atoms has been presented recently [52]. In another attempt, ground and 1s2s 3S, 1S states of a confined He atom were reported [53] taking into account the LDA-approximated exchange-correlation (XC) (with Perdew-Wang parametrization for correlation [54]) with and without self-interaction correction. Response properties such as polarizability and hyperpolarizability of the confined He atom were reported within a DFT-based variation-perturbation approach [55]. In a recent work [56], spherically confined atoms were treated by means of the local exchange potential corresponding to the Zhao-Morrison-Parr and Becke-Johnson potential. Moreover, spherical confinement was used in the comparative study (taking the free-ion limit as reference) of the behavior of spin potential and the pairing energy of first row transition metal cations within the KS model [57]. A detailed analysis of the correlation energy, the performance of several commonly used functionals, electron density, as well as the XC potential in some constrained atoms have been reported [58,59]. The calculation of the static polarizability of confined He and Ne atoms was done through time-dependent DFT in [60].
Recently, there has been a growing interest in information theoretical analysis of diverse models and realistic systems. They have found wide-spread applications in many branches of physics and chemistry, such as thermodynamics, spectroscopy, quantum mechanics. In chemical physics, typically they can provide valuable information regarding the localization-delocalization, diffusion of atomic orbital, periodic properties, spreading of electron density, correlation energy, etc. Entropic uncertainty measures based on these quantities are arguably the most effective quantifiers of uncertainty, as they do not relate to any specific points of the respective Hilbert space. The present work is particularly concerned with Shannon entropy (S) [61,62], which is the arithmetic mean of uncertainty. Interestingly, S like some of the other measures such as Fisher information, Onicescu energy and Rényi entropy are functionals of density and also characterize density. Many articles have been published to analyze these measures in free systems (e.g., for free He, we refer the reader to a recent article [63]), but in confined quantum systems as treated here, analogous studies are quite limited. Two such reports [64,65] in CHA are available so far. A systematic variation of S with respect to r c in r and p spaces has been presented only recently [66,67] for = 0 , as well as non-zero- states. One finds that confinement affects S more profoundly in the stronger regime. Further, S r increases with the rise in r c and at the very low- r c region (≈0.1), CHA displays exactly the opposite trend from a free H atom ( S r declines with the rise in n keeping l fixed). Usually, the effect of perturbation on higher quantum number states is more pronounced. In confined two-electron isoelectronic series (H , He, Li + , Be 2 + ), S has been reported for only ground states [64] by means of BLYP calculations.
Some reports [52,68,69,70,71,72,73] are available on S in penetrable confinement in atoms. For example, it was observed [68] that, in CHA, up to a certain value of r c , S r decreases with r c . However, for small r c and depending on barrier height, S r may also increase. Apart from constant potential, it was probed [73] for confinements imposed by a dielectric continuum and by isotropic harmonic potential. It was also proposed [72] as an indicator to measure the delocalization of electron density. Ground-state atomic S’s, as a function of the width of confining potential, was calculated by employing the correlated Hylleraas-type wave function in both repulsive and attractive finite potentials [70]. Confinement by an inert geometric planar boundary with finite barrier height has been studied [69] within a Thomas-Fermi-Dirac-Weizsäcker-type DFT framework. Some limited works exist on excited states as well, viz. low-lying singly [70] and doubly [71] excited states. It is worth noting that S values in He are available [74] in selected excited states, such as 1 , 3 S e , 1 , 3 P o and 1 , 3 D e .
The objective of this work is to make a thorough systematic analysis of S in a He-like ion placed inside a spherical cage or radius r c . This is done by invoking DFT within a work function-based exchange potential in conjunction with two correlation functionals, viz., a local, parametrized Wigner-type [75] and somewhat involved nonlinear Lee-Yang-Parr (LYP) [76] functional. The pertinent KS differential equation is solved within the Dirichlet boundary condition by means of the GPS method in an accurate and efficient manner. The electron density as well as S r are calculated from the self-consistent orbitals. The p-space orbitals are constructed from respective r-space orbitals via standard Fourier transform, from which the S p ’s are computed. Variation of S r , S p and the total Shannon entropy sum ( S = S r + S p ) with respect to r c is offered for He, Li + and Be 2 + . Apart from the ground state, we also consider singly excited 3S, 3P and 3D states arising out of configurations 1sns (n = 2–4) 3S, 1snp (n = 2–3) 3P and 1snd (n = 3). As is apparent from the preceding discussion, there is a lack of such results in the literature, especially in excited states, and we attempt to provide them. The article is organized as follows. Section 2 outlines the methodology used. Section 3 discusses the results along with a comparison with available references, while Section 4 makes a few concluding remarks.

2. Methodology

Here, we briefly outline the proposed density functional method for ground and excited states of an arbitrary atom centered inside an impenetrable spherical cavity, followed by the GPS method for the calculation of the eigenvalues and eigenenergies of the corresponding KS equation. It may be noted that the present method has been very successfully used for ground and various excited states (such as singly, doubly, triply excited states corresponding to low- and high-lying excitation, valence and core excitation, autoionizing states, hollow and doubly hollow states, very high-lying Rydberg states, satellites states, etc.) of free or unconfined neutral atoms, as well as positive and negative ions in a series of articles [77,78,79,80,81,82,83]. However, it has never been tested for any confinement studies as intended here. Thus, we present an extension of the method for the purpose of confinement effects. Our focus remains on essential portions, omitting the relevant details, which could be found in the above references.
The starting point is the non-relativistic single-particle time-independent KS equation with imposed confinement, which can be conveniently written as (atomic unit employed unless otherwise mentioned),
H ( r ) ϕ i ( r ) = ϵ i ( r ) ϕ i ( r ) ,
where H is the perturbed KS Hamiltonian, written as,
H ( r ) = 1 2 2 + v e f f ( r ) v e f f ( r ) = v n e ( r ) + ρ ( r ) | r r | d r + δ E x c [ ρ ( r ) ] δ ρ ( r ) + v c o n f ( r ) .
In the above, v n e ( r ) and v c o n f ( r ) signify external electron-nuclear attraction and the effective confining potentials, whereas the second and third terms in the right-hand side denote classical Coulomb (Hartree) repulsion and many-body XC potentials, respectively. The desired confinement effect is built into the system by introducing a potential of the following form ( r c refers to the radius of a spherical enclosure),
v c o n f ( r ) = 0 , r r c + , r > r c .
Though DFT has achieved impressive success in explaining the electronic structure and properties of a many-electron system in the ground state in the past four decades, the calculation of excited-state energies and densities has remained a bottleneck. This is mainly due to the absence of a Hohenberg-Kohn theorem parallel to the ground state, as well as the lack of a suitable XC functional valid for a general excited state. In this work, we employed an accurate work function-based exchange potential, which was physically motivated [84,85]. Accordingly, exchange energy is interpreted as the interaction energy between an electron at r and its Fermi-Coulomb hole charge density ρ x ( r , r ) at r , and given by,
E x [ ρ ( r ) ] = 1 2 ρ ( r ) ρ x ( r , r ) | r r | d r d r .
Assuming that a unique local exchange potential v x ( r ) exists for a given state, it can be defined as the work done in bringing an electron to the point r against the electric field generated by its Fermi-Coulomb hole density, leading to the following form,
v x ( r ) = r E x ( r ) d l ,
where the electric field is expressed as,
E x ( r ) = ρ x ( r , r ) ( r r ) | r r | 3 d r .
The Fermi hole can be written in terms of orbitals as,
ρ x ( r , r ) = γ ( r , r ) 2 2 ρ ( r ) ,
where γ ( r , r ) = i ϕ i * ( r ) ϕ i ( r ) is the single particle density matrix and ρ ( r ) is the electron density, expressed in terms of occupied atomic orbitals ( n i denotes occupation number) as,
ρ ( r ) = i = 1 N n i | ϕ i ( r ) | 2 .
While the exchange potential v x ( r ) corresponding to a given state arising from an electronic configuration can be accurately calculated by the above procedure as delineated, the correlation potential v c ( r ) is unknown and must be approximated for practical calculations. The current work incorporates two correlation functionals, namely a Wigner-type [75] and LYP [76]. These two functionals were chosen on the basis of their success in the unconfined atomic excited states, which were recorded in [77,78,79,80,81,82,83]. This work will help shed some light on the applicability of such functionals in the context of confined quantum systems, including those studied here.
With this choice of v x ( r ) and v c ( r ) , the resulting KS differential equation,
1 2 2 + v e f f ( r ) ϕ i ( r ) = ε i ϕ i ( r ) ,
needs to be solved, where v e f f ( r ) is as defined in Equation (2), maintaining the Dirichlet boundary condition. For an accurate and efficient solution, we adopted the GPS scheme leading to a non-uniform, optimal spatial discretization. It is a simple, but very effective method; its success has been demonstrated for many static and dynamic properties of a variety of singular and non-singular potentials of physical and chemical interest [80,81,82,83,86,87,88,89] such as, Coulomb, Húlthen, Yukawa, logarithmic, spiked oscillator, Hellmann potential, etc., along with its recent extension to quantum confinement [15,21,22]. As the method is very well established and documented, in the following, we will mention a very brief summary of it; the details are available in the cited references.
The key characteristic of this approach is to approximate an exact function f ( x ) defined in the interval [ 1 , 1 ] by an N th -order polynomial f N ( x ) ,
f ( x ) f N ( x ) = j = 0 N f ( x j ) g j ( x ) ,
which ensures that the approximation is exact at the collocation points x j ,
f N ( x j ) = f ( x j ) .
Here, we utilize the Legendre pseudo-spectral method where x 0 = 1 , x N = 1 , while x j ( j = 1 , , N 1 ) ’s are defined by the roots of the first derivative of Legendre polynomial P N ( x ) , with respect to x, namely,
P N ( x j ) = 0 .
In Equation (10), g j ( x ) are termed cardinal functions, and as such are expressed as,
g j ( x ) = 1 N ( N + 1 ) P N ( x j ) ( 1 x 2 ) P N ( x ) ( x x j ) ,
fulfilling the unique property that g j ( x j ) = δ j , j . Then, use of a non-linear mapping followed by a symmetrization procedure eventually leads to a symmetric eigenvalue problem, which is solved by standard available software to provide accurate eigenvalues and eigenfunctions.
The p-space wave function is obtained numerically from the Fourier transform of the respective r-space counterpart in the following way,
ξ ( p ) = 1 2 π 3 / 2 ϕ ( r ) e i p . r d r .
It is noted here that ξ ( p ) is not normalized; hence, it needs to be normalized. The normalized r- and p-space densities are represented as ρ ( r ) = i = 1 N n i | ϕ i ( r ) | 2 and Π ( p ) = i = 1 N n i | ξ i ( p ) | 2 , respectively, where n i represents the occupation number of each orbital.
Next S r , S p , and Shannon entropy sum S t are defined in terms of the expectation values of logarithmic probability density functions, which have the forms given below as,
S r = R 3 ρ ( r ) ln [ ρ ( r ) ] d r , S p = R 3 Π ( p ) ln [ Π ( p ) ] d p , S t = S r + S p 3 ( 1 + ln π ) , in 3 dimension .
Here, ρ ( r ) and Π ( p ) are both normalized to unity.
All the computations are done numerically. The convergence is ensured by carrying out calculations with respect to variation in grid parameters, such as the total number of radial points and the maximum range of the grid. It is generally observed that convergence is achieved relatively easily in the lower r c region compared to the r c limit. All results given in the following tables and plots were checked for the above convergence.

3. Result and Discussion

At the onset, it would be appropriate to mention a few points to facilitate the discussion. The net information measures in r and p space of confined many-electron system consist of (i) radial and (ii) angular contributions. The angular part remains invariant in both spaces with respect to change in boundary condition resulting from confinement. Note that, an analogous energy analysis with respect to r c in these and other excited states, would be presented elsewhere. However, it suffices here to mention that, energies for both ground and excited states obtained with the present method are in good agreement with those available in literature. Here our aim is to investigate S for confined two-electron atomic systems at various r c in several low-lying states. These are calculated for spherically confined He and extended to confined iso-electronic members, namely, Li + and Be 2 + . Apart from ground state, following low-lying singly excited states have been considered, viz., 1sns 3S with n= 2–4; 1snp 3P with n = 2–3; 1snd 3D having n = 3. The necessary results are presented in the following tables and plots along with available literature values, for a comparative discussion. All calculations are performed with unit-normalized density.
To start with, Table 1 imprints the numerical results of S r , S p , S t for confined He, Li + and Be 2 + in their ground states. In order to put things in perspective, here and in all other tables, three sets of calculations have been performed at each r c , namely, (i) exchange-only (ii) XC with Wigner correlation (iii) XC with LYP correlation. Throughout the article, these three results are denoted by labels X-only, XC-Wigner and XC-LYP respectively. This can give us an idea how X-only and HF results compare and contrast. Moreover, it will help us in getting a sense of correlation contribution in current context, approximated by two functionals. In all three occasions (i)–(iii), S r ’s increase with rise in r c and finally merge to corresponding free atom entropy at a sufficiently large r c . The reported values of S r in this table further reinforces the previous conclusions [64,72] that impenetrable walls impose confinements in a way that localizes the electron density, and consequently S r when r c 0 . However, it is important to point out that, in all these cases total energy monotonically decreases with r c eventually reaching the free-atom limit. Actually, with reduction in r c the r-space electron density gets compressed and as a consequence, S r decreases. On the contrary, S p gradually abates with progress in r c . At all r c ’s, however, S t maintains the lower bound (6.434) governed by the well-known BBM inequality [62]. The qualitative behavior of S r , S p , S t with growth in r c does not change much with atomic charge (Z), although their numerical values differ. In fact, at a given r c , S r regresses and S p progresses with advancement of Z. With rise in Z, electron density gets compressed and hence such a pattern is noticed. Interestingly, while in one-electron systems S t does not depend on Z, in a many-electron atom, with change of Z it varies [90]. Earlier S t has been mentioned as a measure of correlation in free systems. Our work establishes the same fact in confinement as well. It is noticed that, for all the three species, S’s are identical at very low r c region (≤0.5), without or with (either Wigner or LYP) correlation. Furthermore, these two results begin to differ at larger r c indicating correlation effects to assume more significance in the respective free-atom case. In other words, this implies that, at smaller r c region, XC effect is minimum, which enhances with rise in r c . Similar conclusions have been found in the energy analysis of confined He in [55]. For He and Li + , these have been estimated by BLYP calculation [64] in most of the r c ’s considered here, which are appropriately quoted. Note that the HF values [91] for S r , S p , S t in ground state of unconfined He match reasonably well with our X-only results. Since we are unable to find reference theoretical results for X-only S’s in the hard confinement, for direct comparison, as a matter of check, a couple of comparisons on respective free systems is provided here. Thus the literature S r , S p , S t values of He and Li + within HF method [91,92], employing N-normalized densities, given in footnote of the table, are in reasonable agreement with our X-only values. Recently, in a penetrable confinement calculation within HF, some results on S r have been presented [72]. The same has also been calculated from a DFT-based study with hybrid exchange functional [52]. Results for He in free-limit from both these studies, presented in footnote show quite decent agreement with ours. It may be noted that in these two aforementioned references S r in U 0 corresponds to the impenetrable confinement. Highly accurate benchmark-quality result for S r was calculated from a Hylleraas-type variational method producing a value of 2.7051028 and 1.2552726 for free He and Li + [93] respectively. Correlated results obtained from variational Monte carlo and diffusion Monte Carlo methods [92] for He and Li + are also cited in the footnote. The current single-determinantal approach quite nicely compares with reference values in the table–XC-LYP providing a slight edge over XC-Wigner. It would be worthwhile to make a comparative energy analysis of these two functionals, which we intend to do in near future. No reference entropies are available for Be 2 + .
A careful examination of Table 1 reveals that, X-only, XC-Wigner and XC-LYP results provide similar qualitative trend with respect to changes in r c , for all three species. Hence in Figure 1, X-only S r , S p and S t are plotted for ground state of all three isoelectronic members, as functions of r c in panels (a)–(c). The first two panels imply that, for a fixed Z, S r , S p go up and down respectively with rise in r c . On the contrary, for a given r c , variation of these two quantities with Z shows opposite trend; the former decays and latter develops as Z advances. These results reinforce the inferences drawn from Table 1. Another point to be noted here is that, with lowering in r c the difference between S r and S p corresponding to two successive members of the isoelectronic series, diminishes; in other words, as r c enhances, so does the difference. As r c declines, both average electron-nucleus and electron-electron distances fall down. In stronger confinement regime, the effect of Z on ground state gets dominated by confining potential, resulting in the fact that the three S r , S p plots very nearly coincide. For a fixed Z, variation of S t with r c in panel (c) suggests that the entropy sum reduces dramatically from its free atomic value as r c is lowered. With enhanced confinement, a distinct minimum followed by a maximum shows up in the curve for all He-like ions. The minimum tends to shift towards left as Z progresses. This observed pattern is in consonance with that found in [64].
Now, we move on to some low-lying singly excited state- S r , S p , S t for He, Li + and Be 2 + in Table 2, Table 3 and Table 4. Following the presentation strategy of previous table, it reports results for 1s2s 3S, 1s2p 3P and 1s3d 3D states respectively. Like the ground state, here also in all three states, both X-only and correlation-included S r , S p moves up and down respectively with growth in r c . In all cases again S t obeys the stipulated lower bound [62]. The qualitative pattern of S r , S p , S t with progress in r c remains invariant with change of Z. However, their numerical values alter substantially. Similar to the ground state, at a fixed r c , S r falls off and S p enhances with advancement of Z. Once again, for a given Z, at low r c region, X-only, XC-Wigner and XC-LYP results practically merge with each other, signifying that the effect of correlation is somewhat less impactful in stronger confinement region, as the confining potential leads the contribution in this scenario. Whereas, with rise in r c , the correlation effect prevails, indicating its importance in free conditions. Except the free atom-limit of S r in 1s2s 3S He, no reference results could be found for any of the confined states and we hope the present work would provide useful guideline in future.
In order to gain a better understanding of the effect of confinement on excited states, S r , S p in compressed He have been plotted in panels (a), (b) of Figure 2, for three triplet singly excited states arising from configuration 1sns, corresponding to n =2–4. Since correlation does not affect the results qualitatively, for this purpose, it suffices to consider X-only results. With this in mind, here and in next figure (Figure 3), only X-only entropies are shown. It is obvious from these plots that, for excited states, as in Figure 1, S r gains with rise in r c , while S p declines. Now for a fixed r c , the behavior of S r with n (in this series), shows interesting pattern. For a large enough value of r c , which corresponds to the free-atom limit of He of the state under consideration, S r progresses as n grows. Though it may not be so apparent from the data presented in respective table or plot, as the maximum range of r c presented here is 10 a.u. This can be concluded from the fact that S r for 1s2s, 1s3s and 1s4s triplet states in the free limit are 5.20, 6.53, 7.43 respectively, signifying a progressive delocalization. But this pattern gets dissolved with reduction in r c and crossing between S r for different states occurs. From the inset plots of panel (a) it is noticed that, at r c 7.43 there is a crossing between S r ( 1 s 2 s ) and S r ( 1 s 4 s ) ; another crossing occurs at r c 6.36 between S r ( 1 s 3 s ) and S r ( 1 s 4 s ) . With further gain in confinement strength, at r c 0.33 and 0.34, crossings take place between S r ( 1 s 2 s ) , S r ( 1 s 4 s ) and S r ( 1 s 2 s ) , S r ( 1 s 3 s ) respectively. Hence, one encounters frequent change in the order arrangement on proceeding from free to strong confined regime. In the strong confinement regime ( r c 0.1 ) the following ordering of entropy holds good: S r ( 1 s 2 s ) > S r ( 1 s 3 s ) > S r ( 1 s 4 s ) . Apparently, there exists an interplay between two competing effects, namely, (i) radial confinement (localization) and (ii) accumulation of nodes and humps with growth in n (delocalization). And these two opposing forces control the ordering of S values of these states. Similarly, in p space also, such crossovers prevail at various r c ’s. However, at r c limit and r c = 0.1 , S p shows opposite trend to that of S r .
As a continuation of the earlier discussion, we present in Figure 3, X-only S r , S p and S t as a functions of r c corresponding to 3S, 3P and 3D states resulting from 1s3s, 1s3p and 1s3d configurations of He, in panels (a)–(c). From panel (a), at r c 10 , we obtain the following ordering of entropy: S r ( 3 S ) > S r ( 3 P ) > S r ( 3 D ) , indicating a drop in fluctuation as one passes from 3S to 3P to 3D. Similar to that in Figure 2, here also, multiple crossovers take place at intermediate and lower r c region and eventually settles with the following sequential order S r ( 3 P ) > S r ( 3 S ) > S r ( 3 D ) at ( r c 0.1 ) , representing strong confinement region. Now in conjugate space, at higher r c region 20 , S p displays an exact opposite trend to that of S r in the free limit, which is depicted in panel (b). This is obvious as the more localized a state is in r space, the more diffused it is in p space. Here also, due to crossover between states this pattern gets dissolved at lower r c , leading to an ordering as S p ( 3 S ) > S r ( 3 P ) > S r ( 3 D ) , at r c 0.1 . Panel (c) portrays the response of S t which verifies that the lower bound is maintained throughout entire confinement region.

4. Conclusions

Shannon information entropy (in r and p spaces) was analyzed for confined He iso-electronic series. Ground and excited states were studied via a simple DFT method, by solving the radial KS equation through a generalized Legendre pseudospectral method. Some attempts are known for S of the free He atom, as well as its confinement within a soft, penetrable boundary. However, to the best of our knowledge, this is the first such systematic study of information in a confined two-electron atom within a rigid, impenetrable spherical cage. Apart from the ground state, several low-lying singly excited triplet states of the iso-electronic series were considered. As the X-only entropies were comparable to their HF counterparts in the free atom limit, it was expected that this would also hold in the confined case as well. The effects of electron correlation were probed through two correlation functionals. For the states considered here, the correlation contribution remained rather small in the low r c regime, assuming greater significance as the latter approached the free atom limit. It was observed that the two correlation functionals offered quite comparable results as far as Shannon entropy was concerned. To get more accurate results, it would be necessary to design/employ proper correlation functionals suited for confined systems.
It was seen that S r amplified and S p declined with the rise in r c , in both ground and excited states under consideration. Besides, for a particular confinement strength, as Z grew, the state of a system became more localized with a consequent drop and rise in S r , S p , respectively. For the two families of states arising out of configurations (a) 1sns 3S (n = 2–4) and (b) 1s3s 3S, 1s3p 3P, 1s3d 3D, in the intermediate and lower r c region, the information entropies showed interesting crossovers and finally reached their free atom limit at certain large r c . In all cases, the S t bound was maintained. The emergence of these novel characteristics of S r , S p , and S t makes such information-centric analyses valuable tools for structure and dynamics under a constrained environment. It would be worthwhile to extend the present study to the case of supposedly more realistic penetrable boundaries. Besides, we are also interested in several other information measures like Fisher information, Onicescu energy, complexity, etc., in such systems. Some of these works may be undertaken in the future.

Author Contributions

S.M. contributed in the methodology, program writing, analysis, writing the manuscript and A.K.R. contributed in all these along with editing the manuscript, funding and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

SM is thankful to IISER Kolkata for a Senior Research Fellowship. Financial support from DAE BRNS, Mumbai (Sanction Order 58/14/03/2019-BRNS), is gratefully acknowledged.

Acknowledgments

SM is thankful to IISER Kolkata for a Senior Research Fellowship. Financial support from DAE BRNS, Mumbai (Sanction Order 58/14/03/2019-BRNS), is gratefully acknowledged. Neetik Mukherjee is thanked for the fruitful discussion. We thank the three anonymous referees for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jaskólski, W. Confined many-electron systems. Phys. Rep. 1996, 271, 1–66. [Google Scholar] [CrossRef]
  2. Dolmatov, V.K.; Baltenkov, A.S.; Connerade, J.-P. Structure and photoionization of confined atoms. Radiat. Phys. Chem. 2004, 2004, 417–433. [Google Scholar] [CrossRef]
  3. Sabin, J.; Brändas, E.; Cruz, S. (Eds.) Advances in Quantum Chemistry; Academic Press: New York, NY, USA, 2009; Volume 57, 58. [Google Scholar]
  4. Sen, K.D. (Ed.) Electronic Structure of Quantum Confined Atoms and Molecules; Springer International Publishing: Cham, Switzerland, 2014. [Google Scholar]
  5. Ley-Koo, E. Recent progress in confined atoms and molecules: Superintegrability and symmetry breakings. Rev. Mex. Fís. 2018, 64, 326–363. [Google Scholar] [CrossRef] [Green Version]
  6. Michels, A.; de Boer, J.; Bijl, A. Remarks concerning molecural interaction and their influence on the polarisability. Physica 1937, 1937, 981–994. [Google Scholar] [CrossRef]
  7. Ludeña, E.V. SCF calculations for hydrogen in a spherical box. J. Chem. Phys. 1977, 1977, 468–470. [Google Scholar] [CrossRef]
  8. Marín, J.L.; Cruz, S.A. Enclosed quantum systems: use of the direct variational method. J. Phys. B 1991, 1991, 2899–2907. [Google Scholar] [CrossRef]
  9. Goldman, S.; Joslin, C. Spectroscopic properties of an isotropically compressed hydrogen atom. J. Phys. Chem. 1992, 1992, 6021–6027. [Google Scholar] [CrossRef]
  10. Aquino, N. Accurate energy eigenvalues for enclosed hydrogen atom within spherical impenetrable boxes. Int. J. Quant. Chem. 1995, 1995, 107–115. [Google Scholar] [CrossRef]
  11. Sen, K.D.; Garza, J.; Vargas, R.; Aquino, N. Static dipole polarizability of shell-confined hydrogen atom. Phys. Lett. A 2002, 2002, 299–304. [Google Scholar] [CrossRef]
  12. Laughlin, C.; Burrows, B.L.; Cohen, M. A hydrogen-like atom confined within an impenetrable spherical box. J. Phys. B 2002, 2002, 701–715. [Google Scholar] [CrossRef]
  13. Laughlin, C. On the dipole polarizability of a hydrogen-like atom confined in an impenetrable spherical box. J. Phys. B 2004, 2004, 4085–4099. [Google Scholar] [CrossRef]
  14. Burrows, B.L.; Cohen, M. Exact solutions for perturbed confined hydrogen atoms: Polarizabilities and nuclear shielding factors. Phys. Rev. A 2005, 2005, 032508. [Google Scholar] [CrossRef]
  15. Sen, K.D.; Roy, A.K. Studies on the 3D confined potentials using generalized pseudospectral approach. Phys. Lett. A 2006, 2006, 112–119. [Google Scholar] [CrossRef] [Green Version]
  16. Burrows, B.L.; Cohen, M. Exact solutions for spherically confined hydrogen-like atoms. Int. J. Quant. Chem. 2006, 2006, 478–484. [Google Scholar] [CrossRef]
  17. Aquino, N.; Campoy, G.; Montgomery, H.E., Jr. Highly accurate solutions for the confined hydrogen atom. Int. J. Quant. Chem. 2007, 2007, 1548–1558. [Google Scholar] [CrossRef]
  18. Baye, D.; Sen, K.D. Confined hydrogen atom by the Lagrange-mesh method: Energies, mean radii, and dynamic polarizabilities. Phys. Rev. E 2008, 2008, 026701. [Google Scholar] [CrossRef]
  19. Ciftci, H.; Hall, R.L.; Saad, N. Study of a confined hydrogen-like atom by the asymptotic iteration method. Int. J. Quant. Chem. 2009, 2009, 931–937. [Google Scholar] [CrossRef] [Green Version]
  20. Montgomery, H.E., Jr.; Sen, K.D. Electron density and its derivatives at the nucleus for spherically confined hydrogen atom. Int. J. Quant. Chem. 2009, 2009, 688–692. [Google Scholar] [CrossRef]
  21. Roy, A.K. Spherical confinement of coulombic systems inside an impenetrable box: H atom and the hulthén potential. Int. J. Quant. Chem. 2015, 2015, 937–947. [Google Scholar] [CrossRef] [Green Version]
  22. Roy, A.K. Critical parameters and spherical confinement of H atom in screened Coulomb potential. Int. J. Quant. Chem. 2016, 2016, 953–960. [Google Scholar] [CrossRef] [Green Version]
  23. Ten Seldam, C.A.; De Groot, S.R. On the ground state of a model for compressed helium. Physica 1952, 1952, 891–904. [Google Scholar] [CrossRef]
  24. Ludeña, E.V. SCF Hartree–Fock calculations of ground state wavefunctions of compressed atoms. J. Chem. Phys. 1978, 1978, 1770–1775. [Google Scholar] [CrossRef]
  25. Garza, J.; Hernández-Pérez, J.M.; Ramírez, J.-Z.; Vargas, R. Basis set effects on the Hartree–Fock description of confined many-electron atoms. J. Phys. B 2012, 2012, 015002. [Google Scholar] [CrossRef]
  26. Ludeña, E.V.; Gregori, M. Configuration interaction calculations for two-electron atoms in a spherical box. J. Chem. Phys. 1979, 1979, 2235–2240. [Google Scholar] [CrossRef]
  27. Rivelino, R.; Vianna, J.D.M. A configuration interaction model to investigate many-electron systems in cavities. J. Phys. B 2001, 34, L645–L650. [Google Scholar] [CrossRef]
  28. Joslin, C.; Goldman, S. Quantum Monte Carlo studies of two-electron atoms constrained in spherical boxes. J. Phys. B 1992, 1992, 1965–1975. [Google Scholar] [CrossRef]
  29. Banerjee, A.; Kamal, C.; Chowdhury, A. Calculation of ground-and excited-state energies of confined helium atom. Phys. Lett. A 2006, 2006, 121–125. [Google Scholar] [CrossRef] [Green Version]
  30. Flores-Riveros, A.; Aquino, N.; Montgomery, H.E., Jr. Spherically compressed helium atom described by perturbative and variational methods. Phys. Lett. A 2010, 2010, 1246–1252. [Google Scholar] [CrossRef]
  31. Le Sech, C.; Banerjee, A. A variational approach to the Dirichlet boundary condition: Application to confined H−, He and Li. J. Phys. B 2011, 2011, 105003. [Google Scholar] [CrossRef]
  32. Si, T.-Y.; Bao, C.-G.; Li, B.-W. Energy spectra of the confined atoms obtained by using B-splines. Commun. Theor. Phys. 2001, 2001, 195–200. [Google Scholar]
  33. Montgomery, H.E., Jr.; Aquino, N.; Flores-Riveros, A. The ground state energy of a helium atom under strong confinement. Phys. Lett. A 2010, 2010, 2044–2047. [Google Scholar] [CrossRef]
  34. Aquino, N.; Flores-Riveros, A.; Rivas-Silva, J.F. The compressed helium atom variationally treated via a correlated Hylleraas wave function. Phys. Lett. A 2003, 2003, 326–336. [Google Scholar] [CrossRef]
  35. Flores-Riveros, A.; Rodríguez-Contreras, A. Compression effects in helium-like atoms (Z = 1, …, 5) constrained by hard spherical walls. Phys. Lett. A 2008, 2008, 6175–6182. [Google Scholar] [CrossRef]
  36. Laughlin, C.; Chu, S.I. A highly accurate study of a helium atom under pressure. J. Phys. A 2009, 2009, 265004. [Google Scholar] [CrossRef]
  37. Wilson, C.L.; Montgomery, H.E., Jr.; Sen, K.D.; Thompson, D.C. Electron correlation energy in confined two—electron systems. Phys. Lett. A 2010, 2010, 4415–4419. [Google Scholar] [CrossRef] [Green Version]
  38. Montgomery, H.E., Jr.; Pupyshev, V.I. Confined helium: excited singlet and triplet states. Phys. Lett. A 2013, 2013, 2880–2883. [Google Scholar] [CrossRef]
  39. Bhattacharyya, S.; Saha, J.K.; Mukherjee, P.K.; Mukherjee, T.K. Precise estimation of the energy levels of two—electron atoms under spherical confinement. Phys. Scr. 2013, 2013, 065305. [Google Scholar] [CrossRef]
  40. Montgomery, H.E., Jr.; Pupyshev, V.I. Confined two—electron systems: excited singlet and triplet S states. Theor. Chem. Acc. 2015, 2015, 1598. [Google Scholar] [CrossRef]
  41. Saha, J.; Bhattacharyya, S.; Mukherjee, T.K. Ritz variational calculation for the singly excited states of compressed two—electron atoms. Int. J. Quantum Chem. 2016, 2016, 1802–1813. [Google Scholar] [CrossRef]
  42. Yakar, Y.; Çakir, B.; Özmen, A. Computation of ionization and various excited state energies ofhelium and helium-like quantum dots. Int. J. Quant. Chem. 2011, 2011, 4139–4149. [Google Scholar] [CrossRef]
  43. Doma, S.B.; El-Gammal, F.N. Application of variational Monte Carlo method to the confined helium atom. J. Theor. Appl. Phys. 2012, 2012, 28. [Google Scholar] [CrossRef] [Green Version]
  44. Sarsa, A.; Le Sech, C. Variational Monte Carlo method with dirichlet boundary conditions: Application to the study of confined systems by impenetrable surfaces with different symmetries. J. Chem. Theory Comput. 2011, 2011, 2786–2794. [Google Scholar] [CrossRef] [PubMed]
  45. Young, T.D.; Vargas, R.; Garza, J. A Hartree–Fock study of the confined helium atom: Local and global basis set approaches. Phys. Lett. A 2016, 2016, 712–717. [Google Scholar] [CrossRef]
  46. Pupyshev, V.I.; Montgomery, H.E., Jr. One-and multiconfigurational study of excited states of He atom in a small impenetrable cavity. Theor. Chem. Acc. 2017, 2017, 138. [Google Scholar] [CrossRef]
  47. Parr, R.G.; Yang, W. Density Functional Theory of Atoms and Molecules; Oxford University Press: New York, NY, USA, 1989. [Google Scholar]
  48. Fiolhais, C.; Nogueira, F.; Marques, M. A Primer in Density Functional Theory; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar]
  49. Engel, E.; Dreizler, R.M. Density Functional Theory: An Advance Course (Theoretical and Mathematical Physics); Springer: New York, NY, USA, 2011. [Google Scholar]
  50. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 1988, 3098–3100. [Google Scholar] [CrossRef]
  51. Garza, J.; Vargas, R.; Vela, A. Numerical self-consistent-field method to solve the Kohn-Sham equations in confined many-electron atoms. Phys. Rev. E 1998, 1998, 3949–3954. [Google Scholar] [CrossRef]
  52. Duarte-Alcaráz, F.-A.; Martínez-Sánchez, M.-A.; Rivera-Almazo, M.; Vargas, R.; Rosas-Burgos, R.-A.; Garza, J. Testing one-parameter hybrid exchange functionals in confined atomic systems. J. Phys. B 2019, 2019, 135002. [Google Scholar] [CrossRef]
  53. Aquino, N.; Flores-Riveros, A.; Rivas-Silva, J.F.; Sen, K.D. Confined helium atom low-lying S states analyzed through correlated Hylleraas wave functions and the Kohn–Sham model. J. Chem. Phys. 2006, 2006, 054311. [Google Scholar] [CrossRef]
  54. Perdew, J.P.; Wang, Y. Accurate and simple analytic representation of the electron–gas correlation energy. Phys. Rev. B 1992, 1992, 13244–13249. [Google Scholar] [CrossRef]
  55. Waugh, S.; Chowdhury, A.; Banerjee, A. On the variation of polarizability and hyperpolarizability of a confined atom with the strength of confinement: A case study of a helium atom. J. Phys. B 2010, 2010, 225002. [Google Scholar] [CrossRef]
  56. Vyboishchikov, S.F. Modeling Exact Exchange Potential in Spherically Confined Atoms. J. Comput. Chem. 2015, 2015, 2037–2043. [Google Scholar] [CrossRef] [PubMed]
  57. Lozano-Espinosa, M.; Garza, J.; Galván, M. Confinement effects on the spin potential of first row transition metal cations. Philos. Mag. 2017, 2017, 284–297. [Google Scholar] [CrossRef]
  58. Vyboishchikov, S.F. Correlation energy, correlated electron density, and exchange-correlation potential in some spherically confined atoms. J. Comput. Chem. 2016, 2016, 2677–2686. [Google Scholar] [CrossRef] [PubMed]
  59. Vyboishchikov, S.F. A Simple Local Correlation Energy Functional for Spherically Confined Atoms from ab Initio Correlation Energy Density. ChemPhysChem 2017, 2017, 3478–3484. [Google Scholar] [CrossRef]
  60. van Faassen, M. Atoms in boxes: From confined atoms to electron-atom scattering. J. Chem. Phys. 2009, 2009, 104108. [Google Scholar]
  61. Shannon, C.E. Prediction and Entropy of Printed English. Bell Sys. Tech. J. 1951, 1951, 50–64. [Google Scholar] [CrossRef]
  62. Bialynicki-Birula, I.; Mycielski, J. Uncertainty relations for information entropy in wave mechanics. Commun. Math. Phys. 1975, 1975, 129–132. [Google Scholar] [CrossRef]
  63. Ou, J.-H.; Ho, Y.K. Benchmark calculations of Rényi, Tsallis entropies, and Onicescu information energy for ground state helium using correlated Hylleraas wave functions. Int. J. Quantum Chem. 2019, 119, e25928. [Google Scholar] [CrossRef]
  64. Sen, K.D. Characteristic features of Shannon information entropy of confined atoms. J. Chem. Phys. 2005, 2005, 074110. [Google Scholar] [CrossRef]
  65. Jiao, L.G.; Zan, L.R.; Zhang, Y.Z.; Ho, Y.K. Benchmark values of Shannon entropy for spherically confined hydrogen atom. Int. J. Quant. Chem. 2017, 117, e25375. [Google Scholar] [CrossRef]
  66. Mukherjee, N.; Roy, A.K. Information-entropic measures in free and confined hydrogen atom. Int. J. Quant. Chem. 2018, 118, e25596. [Google Scholar] [CrossRef]
  67. Mukherjee, N.; Roy, A.K. Information-entropic measures for non-zero l states of confined hydrogen-like ions. Eur. Phys. J. D 2018, 2018, 118. [Google Scholar]
  68. Aquino, N.; Flores-Riveros, A.; Rivas-Silva, J.F. Shannon and Fisher entropies for a hydrogen atom under soft spherical confinement. Phys. Lett. A 2013, 2013, 2062–2068. [Google Scholar] [CrossRef]
  69. Cruz, C.; Díaz-García, S.A.; Olivares-Pilón, H.; Cabrera-Trujillo, R. Many-electron atom confinement by a penetrable planar boundary. Radiat. Effects Defects Solids 2016, 2016, 123–134. [Google Scholar] [CrossRef]
  70. Ou, J.-H.; Ho, Y.K. Shannon information entropy in position space for the ground and singly excited states of helium with finite confinements. Atoms 2017, 2017, 15. [Google Scholar] [CrossRef] [Green Version]
  71. Ou, J.-H.; Ho, Y.K. Shannon information entropy in position space for doubly excited states of helium with finite confinements. Chem. Phys. Lett. 2017, 2017, 116–120. [Google Scholar] [CrossRef]
  72. Rodriguez-Bautista, M.; Vargas, R.; Aquino, N.; Garza, J. Electron-density delocalization in many-electron atoms confined by penetrable walls: A Hartree–Fock study. Int. J. Quant. Chem. 2018, 118, e25571. [Google Scholar] [CrossRef]
  73. Martínez-Sánchez, M.-A.; Vargas, R.; Garza, J. Shannon Entropy for the Hydrogen Atom Confined by Four Different Potentials. Quantum Rep. 2019, 2019, 208–218. [Google Scholar] [CrossRef] [Green Version]
  74. Restrepo Cuartas, J.P.; Sanz-Vicario, J.L. nformation and entanglement measures applied to the analysis of complexity in doubly excited states of helium. Phys. Rev. A 2015, 2015, 052301. [Google Scholar] [CrossRef] [Green Version]
  75. Brual, G.; Rothstein, S.M. Rare gas interactions using an improved statistical method. J. Chem. Phys. 1978, 1978, 1177–1183. [Google Scholar] [CrossRef]
  76. Lee, C.; Wang, Y.; Parr, R.G. Development of the Colle–Salvetti correlation–energy formula into a functional of the electron density. Phys. Rev. B 1988, 1988, 785–789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Roy, A.K.; Singh, R.; Deb, B.M. Density-functional calculations for doubly excited states of He, Li+, Be2+ and B3+ (1,3Se, 3Po, 1,3De, 1,3Po, 1Ge,). J. Phys. B 1997, 1997, 4763–4782. [Google Scholar] [CrossRef]
  78. Roy, A.K.; Singh, R.; Deb, B.M. Density functional calculations on triply excited states of lithium isoelectronic sequence. Int. J. Quant. Chem. 1997, 1997, 317–332. [Google Scholar] [CrossRef]
  79. Roy, A.K.; Deb, B.M. Atomic inner-shell transitions: A density functional approach. Phys. Lett. A 1997, 1997, 465–471. [Google Scholar] [CrossRef]
  80. Roy, A.K.; Chu, S.I. Density-functional calculations on singly and doubly excited Rydberg states of many-electron atoms. Phys. Rev. A 2002, 2002, 052508. [Google Scholar] [CrossRef] [Green Version]
  81. Roy, A.K. Studies on the hollow states of atomic lithium using a density functional approach. J. Phys. B 2004, 2004, 4369–4386. [Google Scholar] [CrossRef] [Green Version]
  82. Roy, A.K. Density functional studies on the hollow resonances in the Li-isoelectronic sequence (Z = 4–10). J. Phys. B. 2005, 2005, 1591–1605. [Google Scholar] [CrossRef] [Green Version]
  83. Roy, A.K.; Jalbout, A.F. Ground and excited states of Li−, Be∓ through a density-based approach. Chem. Phys. Lett. 2007, 2007, 355–360. [Google Scholar] [CrossRef] [Green Version]
  84. Sahni, V.; Harbola, M. Quantum-Mechanical interpretation of the local many-body potential of density-functional theory. Int. J. Quant. Chem. Symp. 1990, 1990, 569–584. [Google Scholar] [CrossRef]
  85. Sahni, V.; Li, Y.; Harbola, M. Atomic structure in the Pauli-correlated approximation. Phys. Rev. A 1992, 1992, 1434–1448. [Google Scholar] [CrossRef]
  86. Roy, A.K. Calculation of the spiked harmonic oscillators through a generalized pseudospectral method. Phys. Lett. A 2004, 2004, 231–238. [Google Scholar] [CrossRef] [Green Version]
  87. Roy, A.K. Calculation of the bound states of power-law and logarithmic potentials through a generalized pseudospectral method. J. Phys. G 2004, 2004, 269–278. [Google Scholar] [CrossRef] [Green Version]
  88. Roy, A.K. The generalized pseudospectral approach to the bound states of Húlthen and Yukawa potential. Pramana J. Phys. 2005, 2005, 1–15. [Google Scholar] [CrossRef] [Green Version]
  89. Roy, A.K. Studies on some singular potentials in quantum mechanics. Int. J. Quant. Chem. 2005, 2005, 861–870. [Google Scholar] [CrossRef] [Green Version]
  90. Guevara, N.L.; Sagar, R.P.; Esquivel, R.O. Information uncertainty-type inequalities in atomic systems. J. Chem. Phys. 2003, 2003, 7030–7036. [Google Scholar] [CrossRef]
  91. Gadre, S.R.; Sears, S.B.; Chakravorty, S.J.; Bendale, R.D. Some novel characteristics of atomic information entropies. Phys. Rev. A 1985, 1985, 2602–2606. [Google Scholar] [CrossRef]
  92. Amovilli, C.; Floris, F.M. Shannon entropy in atoms: a test for the assessment of density functionals in Kohn–Sham theory. Computation 2018, 2018, 36. [Google Scholar] [CrossRef] [Green Version]
  93. Lin, C.-H.; Ho, Y.K. Shannon information entropy in position space for two-electron atomic systems. Chem. Phys. Lett. 2015, 2015, 261–264. [Google Scholar] [CrossRef] [Green Version]
  94. Lin, C.-H.; Ho, Y.K. Shannon, Rényi, Tsallis Entropies and Onicescu Information Energy for Low-Lying Singly Excited States of Helium. Atoms 2019, 2019, 70. [Google Scholar]
Figure 1. Variation of S r , S p , S t for He-isoelectronic series (Z = 2–4) with r c in (ac).
Figure 1. Variation of S r , S p , S t for He-isoelectronic series (Z = 2–4) with r c in (ac).
Quantumrep 02 00012 g001
Figure 2. Variation of S r , S p in 1sns 3S (n = 2–4) states of confined He with r c . See the text for details.
Figure 2. Variation of S r , S p in 1sns 3S (n = 2–4) states of confined He with r c . See the text for details.
Quantumrep 02 00012 g002
Figure 3. S r , S p , S t in 1s3s 3S, 1s3p 3P, 1s3d 3D states of confined He, against r c . See the text for details.
Figure 3. S r , S p , S t in 1s3s 3S, 1s3p 3P, 1s3d 3D states of confined He, against r c . See the text for details.
Quantumrep 02 00012 g003
Table 1. S r , S p , S t (a.u.) in the ground states of confined He, Li + , and Be 2 + . See the text for details.
Table 1. S r , S p , S t (a.u.) in the ground states of confined He, Li + , and Be 2 + . See the text for details.
Species r c X-OnlyXC-WignerXC-LYPLiterature 1
S r S p S t S r S p S t S r S p S t S r S p S t
0.1−6.253412.8556.5744−6.253412.8556.5744−6.253412.8556.5744
0.3−3.0049.5806.576−3.0049.5806.576−3.0049.5806.576−2.9889.4886.5
0.5−1.5258.0756.550−1.5258.0756.550−1.5258.0756.550−1.4988.0236.525
10.3896.1186.5070.3876.1196.5060.3886.1186.5060.43266.0786.51
He 2,31.41.225.2736.4931.225.2766.491.225.3276.5471.27395.2346.508
21.974.5476.5171.964.5556.5151.964.5496.5092.00974.5196.528
32.504.0616.5612.484.0806.5602.504.0686.5682.52414.0576.581
42.653.9436.5932.633.9696.5992.643.9526.5922.66513.9456.61
52.683.9216.6082.653.956.602.673.9326.6022.70423.9186.622
62.693.9186.6082.663.956.612.683.936.612.71063.9146.625
72.693.9186.6082.663.956.612.683.936.612.71173.9136.625
0.1−6.266512.8606.5935−6.266512.8606.5935−6.266512.8606.5935
0.3−3.0509.6046.554−3.0509.6046.554−3.0509.6046.554−3.0349.5386.504
0.8−0.3926.8906.498−0.3936.8916.498−0.3926.8906.498−0.3536.8496.496
10.126.3766.4960.1216.3786.4990.1226.3776.4990.16596.3356.501
Li+ 4,521.1355.4316.5661.125.4406.561.135.436.561.1745.46.574
2.51.225.3616.5811.215.3736.581.225.366.581.26185.3316.593
31.245.3466.5861.235.3586.581.245.346.581.28785.3136.601
41.255.3436.5931.235.3556.581.245.346.581.29425.3096.603
71.255.3436.5931.235.3556.581.245.346.58
0.1−6.280112.8666.5859−6.280112.8666.5859−6.280112.8666.5859
0.3−3.1029.6366.534−3.1029.6366.534−3.1029.6366.534
0.5−1.7258.2266.501−1.7258.2266.501−1.7258.2266.501
1−0.2296.7486.519−0.2316.7506.519−0.2296.7486.519
Be2+1.50.1916.3736.5640.1866.3786.5640.1906.3746.564
20.276.3116.5810.266.3176.5770.266.3126.572
2.50.286.3056.5850.276.3116.5810.266.3066.566
30.286.3056.5850.276.3116.5810.266.3066.566
70.286.3056.5850.276.3116.5810.266.3066.566
1 DFT calculation [64], using the BLYP functional. 2 For the free atom, using N-normalized HF density; S r : 4.01 [72], 4.06 [91], 4.0100 [92], 4.0107 [52]; S p : 6.45 [91]; S t : 10.52 [91]. Our exchange-only (X-only) results are 4.00, 6.45, and 10.52 respectively. 3 Correlated S r values in the free He atom: (i) unit-normalized density: 2.7051028 [93], 2.705 [94]; (ii) N-normalized density: 4.40106; (variational Monte Carlo), 4.0256 (diffusion Monte Carlo) [92]. Our exchange-correlation (XC) results are 3.93 (XC-Wigner) and 3.96 (XC-Lee-Yang-Parr (LYP)). 4 HF S r value in the free Li + ion using N-normalized density: 1.1023 [92]. Our X-only result is 1.11. 5 Correlated S r value in the Li + ion: (i) unit-normalized density: 1.2552726 [93]; (ii) N-normalized density: 1.1204 (Variational Monte Carlo), [92]. 1.1143 (diffusion Monte Carlo). Our XC results are: 1.08 (XC-Wigner), 1.10 (XC-LYP).
Table 2. S r , S p , S t (a.u.) in 1s2s 3S states of confined He, Li + , and Be 2 + . See the text for details.
Table 2. S r , S p , S t (a.u.) in 1s2s 3S states of confined He, Li + , and Be 2 + . See the text for details.
Species r c X-OnlyXC-WignerXC-LYP
S r S p S t S r S p S t S r S p S t
0.1−6.217214.1907.9728−6.217214.1907.9728−6.217214.1907.9728
0.5−1.44729.3767.9288−1.44729.3767.9288−1.44729.37567.9288
10.5477.3337.8800.5477.3337.8800.5477.3337.880
22.4175.4097.8262.4155.4107.8252.4175.4097.826
He 143.9563.867.8193.9483.877.8183.9533.867.813
64.603.227.824.593.237.824.593.237.82
6.54.713.117.824.693.137.824.693.137.82
7.54.872.947.814.852.967.814.852.977.82
8.54.992.817.804.962.847.804.962.857.81
105.102.697.795.062.727.785.032.757.78
0.1−6.224614.1917.9664−6.224614.1917.9664−6.224614.1917.9664
0.5−1.49059.3927.9015−1.49059.3927.9015−1.49059.3927.9015
10.4427.4017.8430.4427.4017.8430.4427.4017.843
1.51.4816.3357.8161.4806.3367.8161.4816.3367.817
Li+22.1415.6697.8102.1385.6717.8092.1405.6697.809
32.9104.907.8102.9054.907.8052.9094.9057.814
43.324.497.813.3124.497.8023.314.497.80
73.704.077.773.684.097.773.694.097.78
8.53.724.057.773.704.077.773.704.077.77
103.724.057.773.704.077.773.704.077.77
0.1−6.232014.1917.9590−6.232014.1917.9590−6.232014.1917.9590
0.5−1.53769.4157.8774−1.5379.4157.8774−1.53769.4157.8774
10.3227.4997.8210.3227.5007.8220.3227.4997.821
21.8295.9817.8101.8265.9837.8091.8285.9817.809
Be2+32.425.387.802.1415.397.5312.4185.387.798
42.665.127.782.6515.137.7812.6565.127.776
52.735.037.762.725.057.772.725.047.76
62.755.027.772.735.037.762.745.027.76
202.755.027.772.735.037.762.745.027.76
1 Correlated S r value in the free atom is: 5.239 [94].
Table 3. S r , S p , S t (a.u.) in 1s2p 3P states of confined He, Li+, and Be2+. See the text for details.
Table 3. S r , S p , S t (a.u.) in 1s2p 3P states of confined He, Li+, and Be2+. See the text for details.
Species r c X-OnlyXC-WignerXC-LYP
S r S p S t S r S p S t S r S p S t
0.5−1.3928.5967.204−1.3928.5967.204−1.3928.5967.204
0.8−0.0317.2137.182−0.0317.2137.182−0.0317.2137.182
10.606.5707.1700.606.5707.1700.606.5707.170
33.264.047.303.264.057.313.264.047.30
He 154.093.497.584.083.507.584.083.507.58
64.353.347.694.323.367.684.333.357.68
74.543.237.774.513.257.764.513.257.76
7.64.633.187.814.583.217.794.583.217.79
84.683.157.834.653.187.834.643.187.82
104.863.047.904.823.087.904.783.107.88
0.5−1.4378.6207.183−1.4378.6207.183−1.4378.6207.183
0.8−0.1217.2847.163−0.1217.2847.163−0.1217.2847.163
10.4716.6937.1640.4706.6947.1640.476.697.16
1.51.445.917.351.4385.917.3581.445.917.35
Li+22.0095.297.2992.0055.307.3052.0085.307.308
32.654.857.502.644.867.502.6534.867.513
63.334.477.803.314.497.803.314.487.79
83.394.457.843.364.477.833.374.477.84
93.404.457.853.364.477.833.374.477.84
103.404.457.853.374.477.843.374.477.84
0.5−1.4908.6587.168−1.4908.6587.168−1.4908.6587.168
0.8−0.2357.4027.167−0.2357.4027.167−0.2357.4027.167
1.20.6986.5327.230.6976.5347.2310.6986.5327.23
21.5895.8577.4461.5855.8617.4461.5885.8587.446
Be2+2.51.8925.6707.5621.885.677.551.8905.7017.591
52.375.417.782.355.437.782.365.427.78
62.385.417.792.365.427.782.375.427.79
102.385.417.792.365.427.782.375.427.79
1 Correlated S r value in the free atom is: 5.356 [74].
Table 4. S r , S p , S t (a.u.) in 1s3d 3D states of confined He, Li+, and Be2+. See the text for details.
Table 4. S r , S p , S t (a.u.) in 1s3d 3D states of confined He, Li+, and Be2+. See the text for details.
Species r c X-onlyXC-WignerXC-LYP
S r S p S t S r S p S t S r S p S t
0.5−1.30539.1107.8047−1.30539.1107.8047−1.30539.1107.8047
10.7077.0737.7800.7067.0737.7790.7077.0737.780
1.51.8285.9387.7661.8265.9397.7651.8275.9387.765
2.63.1554.6407.7953.1504.6457.7953.1534.6417.794
He 143.9603.9577.9173.9513.9677.9183.9573.9607.917
54.3273.6878.0144.3173.6998.0164.3233.6918.014
64.613.488.094.603.508.104.603.498.09
74.853.338.184.833.348.174.833.348.17
85.043.198.235.033.218.245.023.228.24
105.362.978.335.352.998.345.273.048.24
0.5−1.3399.1297.790−1.3399.1297.790−1.3399.1297.790
0.8−0.0037.7747.771−0.0037.7747.771−0.0037.7747.771
10.6017.1617.7620.6017.1617.7620.6017.1617.762
2.52.6115.277.8812.6085.2807.8882.6115.2777.888
Li+32.9185.0427.9602.9145.0477.9612.9175.0437.960
43.3734.728.0933.364.738.093.374.728.09
53.7044.508.2043.694.518.203.704.508.20
6.54.064.268.324.054.278.324.054.278.32
74.164.208.364.154.218.364.144.218.35
7.54.244.148.384.234.158.384.224.168.38
104.533.978.504.513.998.504.464.038.49
0.5−1.38039.1577.776−1.38039.1577.776−1.38039.1577.776
0.8−0.0977.8577.760−0.0977.8577.760−0.0977.8577.760
10.4607.3037.7630.4597.3047.7630.4607.3037.763
2.52.215.798.002.2115.798.0012.2135.798.003
Be2+32.4965.598.0862.4925.608.0922.4955.598.085
42.9115.328.2312.9075.328.2272.9095.328.229
73.534.948.473.514.968.473.514.968.47
83.604.918.513.594.928.513.584.938.51
8.53.634.908.533.614.928.533.604.928.52
103.664.908.623.644.918.553.624.928.54
1 Correlated S r value in the free atom is: 6.634 [74].

Share and Cite

MDPI and ACS Style

Majumdar, S.; Roy, A.K. Shannon Entropy in Confined He-Like Ions within a Density Functional Formalism. Quantum Rep. 2020, 2, 189-207. https://doi.org/10.3390/quantum2010012

AMA Style

Majumdar S, Roy AK. Shannon Entropy in Confined He-Like Ions within a Density Functional Formalism. Quantum Reports. 2020; 2(1):189-207. https://doi.org/10.3390/quantum2010012

Chicago/Turabian Style

Majumdar, Sangita, and Amlan K. Roy. 2020. "Shannon Entropy in Confined He-Like Ions within a Density Functional Formalism" Quantum Reports 2, no. 1: 189-207. https://doi.org/10.3390/quantum2010012

Article Metrics

Back to TopTop